Decoherence in the dynamical quantum phase transition of the transverse Ising chain 



Sarah Mostame, Gemot Schaller, and Ralf Schiitzhold* 
Institut fur Theoretische Physik, Technische Universitdt Dresden, 01062 Dresden, Germany 

For the prototypical example of the Ising chain in a transverse field, we study the impact of 
decoherence on the sweep through a second-order quantum phase transition. Apart from the advance 
in the general understanding of the dynamics of quantum phase transitions, these findings are 
relevant for adiabatic quantum algorithms due to the similarities between them. It turns out that 
(in contrast to first-order transitions studied previously) the impact of decoherence caused by a weak 
coupling to a rather general environment increases with system size (i.e., number of spins/qubits), 
which might limit the scalability of the system. 
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Recently, the dynamics of quantum phase transitions 
[H attracted increasing interest, see, e.g., @, H, 0]. In. 
contrast to thermal transitions (usually driven by the 
competition between energy and entropy), they are char- 
acterized by a fundamental change of the ground state 
structure (e.g., from para- to ferro-magnetic) at the crit- 
ical value of a variable external parameter (e.g., mag- 
netic field). Quantum phase transitions are induced by 
quantum rather than thermal fluctuations and thus may 
occur at zero temperature. At the critical point, the 
energy levels become arbitrarily close and thus the re- 
sponse times diverge (in the continuum limit). Conse- 
quently, during the sweep trough such a phase transition 
by means of a time-dependent external parameter, small 
external perturbations or internal fluctuations become 
strongly amplified - leading to many interesting effects, 
see, e.g., 0, HI, 0- One of them is the anomalously 
high susceptibility to decoherence (see also Q): Due to 
the convergence of the energy levels at the critical point, 
even low-energy modes of the environment may cause ex- 
citations and thus perturb the system. Here, we study 
the decoherence caused by a small coupling to a rather 
general reservoir for the quantum Ising chain in a trans- 
verse field, which is considered a prototypical example [l| 
for a second-order quantum phase transition (and further 
possesses the advantage of being analytically solvable). 

Apart from the general understanding of quantum 
phase transitions, these investigations are also relevant 
for quantum computing: By constructing the Hamilto- 
nian appropriately, it is possible to encode the solution 
to hard computational problems (such a factoring large 
numbers) in its ground state. In order to reach this solu- 
tion state, we may start off with a simpler Hamiltonian 
whose ground state is easy to prepare as the initial con- 
figuration. If we now steadily transform it to the problem 
Hamiltonian, the adiabatic theorem tells us that we stay 
near the ground state if the evolution is slow enough - 
and thus finally end up in (or close to) the desired solu- 
tion state (adiabatic quantum computing (Icl.[ll]|). How- 
ever, somewhere on the way from the simple initial config- 
uration to the final state, there is typically a critical point 
which bears strong similarities to a quantum phase tran- 



sition (e.g., vanishing gap and diverging entanglement 
in the continuum limit II, 13|). Based on this similar- 
ity, it seems fl3| that adiabatic quantum algorithms cor- 
responding to second-order quantum phase transitions 
should be advantageous compared to isolated avoided 
level crossings (which are analogous to first-order tran- 
sitions). For an adiabatic quantum algorithm (Grover's 
search routin e 11411 ) based on a single isolated avoided 
level crossing j 1 51 ] , the impact of decoherence induced by 
a low-tcmpcrature bath with a well-behaved spectral dis- 
tribution does not destroy the scalability of the system. 
However, as we shall see below, the situation may be 
very different for second-order transitions. These inves- 
tigations are particularly relevant in view of the recent 
announcement (see, e.g., regarding the construction 
of an adiabatic quantum computer with 16 qubits in the 
form of a two-dimensional Ising model. 

The open system under consideration is described by 
the total Hamiltonian H which can be split up into that 
of the closed system H sys and the bath -ffbath acting on 
independent Hilbert spaces Sj sys ® iObath = £> 



H(t) = H sys (t) + ff ba th(i) + A ff in t (t) . 



(1) 



plus an interaction A H mi between the two, which is sup- 
posed to be weak A <C 1 in the sense that it does not 
perturb the state of the system drastically. Note, how- 
ever, that the change of the bath caused by the inter- 
action with the system need not be small. In order to 
describe the evolution of the combined quantum state 
|$(i)) € Sj, we expand it into the instantaneous system 
energy eigenbasis H sys (t)\^ s (t)) = E s (t)\^ s (t)) in 9j sys 
via |$(t)) = J2 S °a(*)l*s(*)) ® l«s(*))) where a s are the 
corresponding amplitudes and \a s ) € .flbath denote the 
associated (normalized but not necessary orthogonal) 
states of the reservoir. Insertion of this expansion into the 
Schrodinger equation i|$(t)) = H{t)\&(t)) yields (h = 1) 







(^ s |g sys i^ r ) 

AE sr 



-i\(a s \(^s\H in t\^r)\ar 



(2) 
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with the energy gaps AE sr (t) = E s (t) — E r (t) of the sys- 
tem and the total phase (including the Berry phase) 

Vs (t) = J* dt'{E s (t>) + H&ktf) + A2T£(f) 

-i^ s (t')\^ s (t'))-i{a s (t')\a s {t'))}, (3) 

containing the energy shift H?^ t = {a s \{^ s \H- lnt \^ r )\a r ) 
and.ff^ th = (a s |i?bath|ov)- Evidently, there are two con- 
tributions for transitions in the Hilbert space fj sys of the 
system: The first term on the right-hand side of Eq. ((2J) 
describes the transitions caused by a non-adiabatic evo- 
lution (see, e.g., [lH). Note, however, that the factor 
(a s \a r ) and the additional phases in Eq. ([3|) give rise to 
extra terms in the adiabatic expansion. The second term 
in Eq. flS)) directly corresponds to transitions caused by 
the interaction with the bath. Since we are mainly inter- 
ested in the impact of the coupling to the bath, we shall 
assume a perfectly adiabatic evolution 17[ of the system 
itself (* s |ij sys |* r ) < AE* r such that the first term in 
Eq. ([2]) is negligible and the second one dominates. 

The quantum Ising chain of n spins we are going to 
study exhibits a time-dependent nearest-neighbor inter- 
action g{t) plus transverse field B(t) = 1 — g(t) 



H sys (t) = - J2 { [1 - <?(*)] + 9(t) crXi+i } 



(4) 



where crj = (crj, crj, o|) are the spin-1/2 Pauli matrices 
acting on the jth qubit and periodic boundary conditions 
<r n+ i = (j\ are imposed. Choosing g(0) = and g(T) = 1 
where T is the evolution time, the system evolves from 
the paramagnetic state | — >— >— > . . . ) trough a second- 
order p| quantum phase transition at g CT = 1/2 to the 
ferromagnetic phase | f"TT ■ • • ) or I III ■ • • )• 

A major advantage of the above Hamiltonian is 
that it can be diagonalized exactly Let us briefly 
review the main steps of the diagonalization of H sys , 
where we switch temporarily to the Heisenberg pic- 
ture for convenience: The set of n qubits in (TjQ) can 
be mapped to a system of n spinless fermions Cj 
via the Jordan- Wigner [2()j transformation given by 
0-7 = 1- 2c]cj and cr? = -(cj + Cj) Y\ e<J o x t . In terms 
of the fermionic operators cj and Cj , the Hamiltonian as- 
sumes a bilinear form containing (1 — g)(l — 2ctcj) and 
g(cj + iCj + cj. +1 Cj + cjcj+i + cjcj +1 ), i.e., the fermion 

number rij — ejej is not conserved in general. This 
bilinear form can now be diagonalized by a Fourier trans- 
formation Cj = ~^2 k Cfc e~ %k ti a ) / i/n followed by a Bogoli- 
ubov 21] transformation c k (t) — Uk(t) ^ k + iv%{t) 7l fe . 
Since the new set of fermionic operators j k is supposed 
to be time-independent, the Bogoliubov coefficients 
u k and v k must satisfy 0] the equations of motion 
iiik = -a k u k + (3 k v k and iv k = a k v k + (3 k u k , where 



a k (t) = 2- 4g(t) cos 2 (fca/2) and (3 k {t) = 2g(t) sin(fca). 
For an adiabatic evolution (^^6^1^,.) <C AE^ r , 
these equations of motion can be solved approximately 
via u fc (t) » [a k (t) + e fc (t)]exp{-i f* dt'e k (t')}/jV k (t) 
as well as v k {t) -f3 k (t)exp{-i dt'e k (t')}/Af k (t) 
with the normalization A4 = yj2e^~\~2a k t k ensuring 
\u k \ 



+ \ v k\ = 1 and the single-particle energies 
e k (t) = 2y/l -Ag(t) [1 - g(t)] cos 2 (fca/2) . 



(5) 



All the excitation energies e& assume their minimum val- 
ues e™ ln = 2|sin(fca/2)| at the critical point g cr = 1/2. 
In the following, we study the scaling of the involved 
quantities in the continuum limit u | oo. In view of 
the fc-spectrum k E tt(1 + 2Z)/(an) : \ka\ < it, 
where a is the lattice spacing, the minimum gap scales 
as AE m [ n — 0(l/n). Finally, the Hamiltonian (H|) reads 



t 1 

7fe7fc - 2 



(0) 



and hence its (instantaneous) ground state contains no 
fermionic quasi-particles y k Jk\^o{t)) — 0. Without the 
environment, the number of fermionic quasi-particles 
l\lk would be conserved and the system would stay in an 
eigenstate (e.g., ground state) for an adiabatic evolution. 

Of course, the impact of decoherence depends on the 
properties of the bath and its interaction with the sys- 
tem (decoherence channels). In order to derive generally 
applicable results, we do not specify the bath -ffbath in 
much detail and start with an interaction A H int which 
is always present: In the Hamiltonian H sys in Eq. (fj[|), 
the transverse field B(t) = 1 — g(t) appears as a clas- 
sical control parameter B c \. However, the external field 
B — » B c \ + SB does also possess (quantum) fluctuations 
SB, which couple to the system of Ising spins. Therefore, 
we start with the following interaction Hamiltonian 



(7) 



where SB denotes the reservoir operator. Incidentally, 
this interaction Hamiltonian yields the same matrix el- 
ements as the non-adiabatic corrections (4 , s |ij sys |^' r ) in 
Eq. ([2]), which can therefore be calculated analogously. 

Starting in the system's ground state ao(t = 0) = 1, 
the excitations s > caused by the weak interaction 
A-ffj n t with the bath 2l s = a s (T) exp{itp s (T)} can be cal- 
culated via response theory, i.e., the solution of Eq. ([2]) 
to first order in A <C 1 is 

a, « -ix [ du)f s {u) j dt(ij s (t)\J2^\Mt)) x 



x exp < i 



-u>t+ f dt' AE s0 (t') 
Jo 



(8) 
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We have subsumed all relevant properties of the environ- 
ment into the spectral function f(uj) of the bath 



dioe-^f^to), (9) 



where A<p s coincides with ip s — ipo in Eq. ([3]) apart from 
the system's energy gap AE s0 and is typically dominated 
by the contribution from H£l th — A^ath- As a ^ rs ^ a P~ 
proximation, we assume that f(u>) does not change sig- 
nificantly if we increase the system size n (scaling limit). 

After inserting the Jordan- Wigner [2(| transformation, 
the matrix element in Eq. ([8]) reads 



£>.|*?(*)ltfo> 



where the w sign refers to the adiabatic approximation. 
Thus, it is only non- vanishing for excited states \ip a ) con- 
taining two quasi-particles s — (k, —k) with opposite mo- 
menta and hence we get AE s0 = 2e^. In order to solve 
the remaining time integrals, it is useful to distinguish 
different w-regimes: First of all, in order to have a quan- 
tum phase transition (or a working adiabatic quantum 
computer), the environment should be cold enough to 
permit the preparation of the system in the initial ground 
state, i.e., uj -C 2 = euit = 0). For intermediate positive 
frequencies 2 ^> uj ^S> A_E™ n ps 2|fca|, we may solve the 
time integral via the saddle-point (or stationary phase) 
approximation. The saddle-point condition for the expo- 
nent in Eq. © reads uj — AE s0 (t* ) = 2tk (t r ) , which cor- 
responds to energy conservation. This condition yields 
two saddle points shortly before and after the transition 
g(tf) ps 1/2 ± \/uj 2 - 4fc 2 a 2 /8. For the spectral excita- 
tion amplitude 21" defined via 2l s — J doj f s (ui) 21" , the 
saddle-point approximation yields 



21 



= o 



Xka 



(11) 



Of course, the result depends on the interpolation dy- 
namics git), see Table UJ For a constant speed interpo- 
lation g(t) = t/T, the run-time needed for an adiabatic 
evolution scales as T — 0(n 2 ) due to the minimum gap 
AE m i n = 0(l/n). For adapted interpolation dynam- 
ics g(t) tx AE(t) or g{t) tx AE 2 it), however, one may 
achieve shorter run-times of T = O(nlnn) or T — 0(n), 
respectively [lit ]. 

The next (higher-order) terms of the saddle-point ex- 
pansion scale with 0{\g{t*)[uj 2 - 4fc 2 a 2 ]~ 1 / 2 /cj) and 
hence the saddle-point approximation breaks down if u) 
approaches the minimum gap AE™ n w 2|fca|, see Fig.Q] 
In this case, we may obtain an upper bound for the 
time integral in Eq. (JSj> via omitting all phases, see Ta- 
ble U For frequencies far below 2|fca|, the saddle points 
at g{t^) ~ 1/2 ± \]u> 2 — 4fc 2 a 2 /8 move away from the 





1 > ui > 2ka 


1 > uj « 2ka 


m = o 


0(\kauj- 1 n) 


0{\n 2 Lu\nu) 


g(t) oc AE(t) 


0(\kauj- 3/2 ^n) 


0(Xn Inn) 


g(t)xAE 2 (t) 


0(Xkauj' 2 ) 


G{Xn) 



TABLE I: Scaling of the spectral excitation amplitude 21" 
in the saddle-point approximation (uj 3> 2ka) and its upper 
bound (uj ~ 2ka) for different interpolation dynamics g(t), 
where AE(t) — 2ej_ 7r /( an ) (t) denotes the fundamental gap. 
In all cases, the total excitation probability (sum over all ui 
and k) increases with system size n. 
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FIG. 1: [Color online] Sketch of the excitation spectrum of 
the Ising chain H aya as a function of g. For a given frequency 
uj > 0, real saddle points correspond to intersections of the 
(solid) energy level curves (e.g., AEi) with the (dashed) ver- 
tical oj-line which occur shortly before (g*) and after (gt) 
the quantum phase transition at g CT = 1/2. The saddle-point 
approximation can only be applied if the intersection angle is 
large enough, i.e., for the drawn uj > line, it would work for 
AEi, but not for AE 3 etc. 



real axis and thus the amplitudes are exponentially sup- 
pressed in the adiabatic limit 



21 



w<C2| ka\ 



D(Xexp{-T(ka) 2 /2}) 



(12) 



for g(i) — t/T and similarly for the other interpolations. 
Finally, for negative frequencies uj < 0, the saddle points 
collide with the branch cut generated by the square- 
root in efe. By deforming the integration contour into 
the complex plane.it can be shown via an argumenta- 
tion analogous to [l8| that the amplitudes are also ex- 
ponentially suppressed in this case. This result can be 
understood in the following way: For frequencies uj be- 
low the lowest excitation energies, the energy uj of the 
reservoir modes is not sufficient for exciting the system 
via energy-conserving transitions. Hence excitations can 
only occur via non-adiabatic processes for which energy- 
conservation becomes ill-defined, but these processes are 
suppressed if the evolution is slow enough. 

In summary, we studied the quantum phase transition 
from paramagnetic to ferromagnetic phase in the quan- 
tum Ising chain in a transverse field via its analytical 



4 



diagonalization and calculated the excitation probabili- 
ties [17| caused by a weak coupling to a rather general 
environment (including possible non-perturbative behav- 
ior of the reservoir). Since the Ising model is considered 
[l| a prototypical example for a second-order quantum 
phase transition, we expect our results to reflect general 
features of second-order transitions. For the decoherence 
channel which is always present (though possibly not 
the dominant channel), we already found that the to- 
tal excitation probability increases with system size n 
(continuum limit): Even though the probability for the 
lowest excitation k — ±7r/(an) can be kept under con- 
trol for a bath which is well-behaved in the infra-red 
limit (see also |15|). the existence of many excited states 
k G 7r(l + 2Z)/(cm) : \ka\ < ir converging near the 
critical point causes the growth of the error probability 
for large systems. This growth can be slowed down a bit 
via adapted interpolation schemes g(t), but not stopped. 
Other decoherence channels will display the same general 
behavior: E.g., for w ^> \ka\, the associated amplitudes 
scale as 21" = ©(A^s^,)/-^^*)), where 4> s denotes the 
matrix element in analogy to (|10p . Typically, for a ho- 
mogeneous coupling to the bath, 4> s does not strongly 
depend on the system size n (for given ka and oj). Since 
<?(i*) decreases for n | oo or at least remains constant [for 
g(t) oc AE 2 (t)], the total excitation probability again in- 
creases with system size n [23j |. 

Using the analogy between adiabatic quantum algo- 
rithms and quantum phase transitions 1^, ll| , this result 
suggests scalability problems of the corresponding adia- 
batic quantum algorithm - unless the temperature of the 
bath stays below the (n-dependent) minimum gap (lij 
or the coupling to the bath decreases with increasing n. 
These problems are caused by the accumulation of many 
levels at the critical point g = 1/2, which presents the 
main difference to isolated avoided level crossings (corre- 
sponding to first-order phase transitions) discussed ear- 
lier [15J. It also causes some difficulties for the idea of 
thermally assisted quantum computation (see, e.g., [22j]) 
since, in the presence of too many available levels, the 
probability of hitting the ground state becomes small. 

Therefore, in order to construct a scalable adiabatic 
quantum algorithm in analogy to the Ising model, suit- 
able error-correction methods will be required. As one 
possibility, one might exploit the quantum Zeno effect 
and suppress transitions in the system by constantly 
measuring the energy, see for example [24j . As another 
interesting idea, let us study a spatial sweep through 
the phase transition, i.e., we do not cross the criti- 
cal point in a homogeneous way, but adopt the fol- 
lowing step-wise interpolation: Starting from the ini- 
tial Hamiltonian af + erf + °f + °"i + • ■ ■ j we change 
it slowly to <j\<j\ + erf + erf + . . . and afterwards to 
a i a 2 + cr 2 <T 3 + a 4 + • • ■ etc. This corresponds to a nonlin- 
ear interpolation path between the two Hamiltonians. In 
this case, the minimum gap (in the relevant subspace that 



is even under bit flip) remains independent of the system 
size n and the run-time T scales linear in n (number of 
steps). Hence, decoherence could be strongly suppressed 
for a low-temperature bath. Of course, the generalization 
of all these concepts and results to more interesting cases 
such as the (NP-complete) two-dimensional Ising model 
is highly non-trivial and requires further investigations. 
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